The Breakdown of Linear Elastic Fracture Mechanics near the Tip of a Rapid Crack 
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We present high resolution measurements of the displacement and strain fields near the tip of a 
dynamic (Mode I) crack. The experiments are performed on polyacrylamide gels, brittle elastomers 
whose fracture dynamics mirror those of typical brittle amorphous materials. Over a wide range 
of propagation velocities (0.2 — 0.8cs), we compare linear elastic fracture mechanics (LEFM) to the 
measured near-tip fields. We find that, sufficiently near the tip, the measured stress intensity factor 
appears to be non-unique, the crack tip significantly deviates from its predicted parabolic form, and 
the strains ahead of the tip are more singular than the r~^/^ divergence predicted by LEFM. These 
results show how LEFM breaks down as the crack tip is approached. 

PACS numbers: 46.50.+a, 62.20.Mk, 89.75.Kd 



The behavior of a cracked body under applied stress is 
of extreme practical and fundamental importance. The 
accepted approach to describing crack dynamics is linear 
elastic fracture mechanics (LEFM) [1]. This framework 
assumes that a material is linearly elastic with all non- 
linear and dissipative processes well confined to the near 
tip vicinity. LEFM provides a full description of the elas- 
tic fields surrounding the tip of a single crack, whether 
static or propagating. Open questions, however, such 
as criteria for crack path selection [2,], the origin of dy- 
namic instabilities (micro-branching) [2| and oscillations 
of a single crack at high velocities [3] underline the need 
for a more detailed understanding of the near tip behav- 
ior in dynamic fracture. Studies of this region have been 
limited mainly to static cracks [J, |5[ , as detailed measure- 
ment of a microscopic region which may move at speeds 
approaching sound speeds entails enormous difficulties. 

Here we report on recent experimental results in the 
dynamic fracture of brittle polyacrylamide gels. In these 
soft materials, fracture dynamics are identical to those 
observed in standard brittle amorphous materials, but 
crack velocities, which scale with material sound speeds, 
are nearly three orders of magnitude lower [6[. Thus, 
detailed visualization of the fields in the near-tip region 
within an effectively 2D medium becomes possible. The 
measurements provide a detailed description of the near- 
tip fields. Quantitative comparison with LEFM delin- 
eates its domain of validity. These findings suggest the 
importance of nonlinear elasticity as the crack tip is ap- 
proached. 

The brittle gels used in these experiments are transpar- 
ent, neo-Hookean, incompressible elastomers that were 
used in [3[. They are composed of 13.8% total monomer 
and 2.6% bis-acrylamide cross-linker concentrations. The 
shear (/i = 35.2±1.4kPa) and Young's {E = 3fi) moduli 
of these gels yield shear and longitudinal wave speeds of 
C5 = 5.90±0.15m/s and cl = 11.8±0.3m/s. Their typical 
dimensions were {XxYxZ) 130 x 125 x 0.2mm where X, 
F, and Z are, respectively, the propagation, loading, and 
gel thickness directions. Reducing the gel thickness sup- 



presses micro-branching and enables single-crack states 
to attain high velocities in effectively 2D media [3| . 

The gels were cast between two flat glass plates. The 
face of one plate was randomly scratched with No. 600 
Alundum lapping powder. These scratches, of 16/im 
mean depth, were imprinted on one of the gel faces. The 
resulting scratch pattern was used as a "tracer field" for 
visualization of the displacement field, and did not affect 
the crack dynamics. 

Experiments were performed as in [6] by imposing 
uniaxial (Mode I) loading via constant displacement in 
the vertical (F) direction. Once the desired stress was 
reached, a seed crack was imposed at the sample's edge, 
midway between the vertical boundaries. As Fig. 1(a) 
shows, the displacement field of the dynamic crack that 
emerged was visualized at the center of the gel by a high- 
speed (485Hz frame rate) camera that was focused on a 
fixed 15.4mm (X) by 12.3mm (Y) region. The camera's 
1280x1024 pixel resolution enabled a 12/im spatial reso- 
lution. Image blur was limited to less than a single pixel 
by using a 2.8/is exposure time. 

A fundamental quantity of fracture mechanics is the 
displacement field n, defined by x^ = x-\-u{x)^ where 
X is an un-deformed 'rest' configuration and x' a de- 
formed one. We measured this quantity as follows. We 
first visualized the scratch pattern by shadowgraphy [7( . 
We then found the displacement field between a pho- 
tograph taken immediately preceding fracture initiation 
(reference frame), and photographs capturing traveling 
cracks. This was done by means of a particle track- 
ing technique where the scratches were tracked by cross- 
correlating small boxes (10-20 pixels in length, at 5 pixel 
intervals) from the reference frame with corresponding 
regions deformed by the crack's passage. The maximum 
correlation for each box provided a sub-pixel measure- 
ment of the displacement field generated by the crack. 
Correcting for the uniform stretch of the reference frame 
yielded u{x) in the rest frame, thereby enabling a direct 
comparison to LEFM. 

Let us first consider the crack tip opening displacement 
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FIG. 1: (Color online) (a) (left) A single crack traveling from 
left to right in a (to scale) uniaxially loaded sample, (center) 
Photograph (here 5.4 x 5.2mm) of the profile at the tip of a 
dynamic crack (v^0.73cs). A parabolic fit (dashed red curve) 
agrees well with the crack profile for a distance larger than 
5mm. (right) A closeup view (1.0x1.8mm) reveals a deviation, 
(5, between the tip of the parabola and the crack tip. (b) S 
and (c) the parabola curvature, a, as functions of v/cg. 



(CTOD), which is the clearest manifestation of the dis- 
placement field in the cracked body. LEFM predicts a 
parabolic CTOD: 
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We define moving frame coordinates, (r, 0) (r = 
y^{x — vty + y'^ and = id.ii~^[y/ {x — vt)])^ centered at 
the crack tip where ^ = is the propagation direction. 
In Eq.([T]), Kj is the stress intensity factor, T is the sub- 
leading correction known as the "T-stress" and Qx,y{^'-, '^) 
are universal functions [8] of and the crack velocity, v. 
As T/E <C 1, the curvature, a, of the parabola provides 
a direct measurement of Kj. This important quantity, 
according to LEFM, wholly dictates the behavior of a 
moving crack. 

To obtain i^/, we extracted the crack opening pro- 
files from the photographs and fitted them to parabolic 
forms over l-2mm from the crack tip, which is within 
the range where LEFM might be expected to be rele- 
vant. The parabolic profiles obtained from the fitted data 
agreed well with CTOD profiles at significantly larger 
scales, indicating their robustness. For example, the fit- 
ted parabola in Fig. 1(a) is indistinguishable from the 
CTOD profile over 5mm. The crack-tip curvature is seen 
to be a decreasing function of the velocity [Fig. 1(c)]. In- 
terestingly, it is well described by a linear function with 
its extrapolated zero intercept at a velocity close to Cg. 

A closer look, however, [Fig. 1(a) right] reveals that a 
noticeable deviation between the fit and the crack edges 
occurs in the near vicinity of the tip. This deviation, 
(5, increases rapidly with crack velocity, from ~ 30/im at 
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FIG. 2: (Color online) Scratch patterns (inset) imprinted on 
the gels are used to measure the displacement field surround- 
ing the crack tip via particle tracking. The measured Ux 
(v = 0.37cs) along the ^ = axis of symmetry (dashed white 
line) is fitted to Eq. (|2]). The values of Ux are accurate up to 
a small additive constant (< 0.05mm). 



low velocities (~ 0.2cs) up to ~ 300/im at high veloci- 
ties (~ 0.8cs) [Fig. 1(b)]. The near-tip region is one of 
considerably higher curvature than the "far-field" regions 
captured by the parabolic form. Attempts to force the tip 
of the parabola to coincide with the measured crack tip 
yielded terrible compatibility with the data at all scales. 
We are therefore led to the conclusion that the near-tip 
region described by S is not described by LEFM. We will 
use 5 to characterize the scale of this divergence. Al- 
though somewhat arbitrary, other length-scales (e.g. the 
distance between the crack tip and the point where the 
divergence begins) display similar functional behavior. 

The near-tip region characterized by 5 can not be ex- 
plained in the framework of LEFM. A natural explana- 
tion for this region of deviation would be to associate 
it with the "process zone" scale, where elasticity breaks 
down due to extreme stresses, giving way to plastic defor- 
mations and fracture itself. LEFM avoids any treatment 
of this zone and regards it as a singular point under the 
conditions of small-scale yielding [l|. Another possible 
signature of the process zone may be the white region 
visible at the tip of the crack [Fig. 1(a) right]. In this re- 
gion, the high strains induce large material deformations 
giving rise to lensing effects. In pictures, as this one, 
where the camera's focal plane is set slightly below the 
gel (between the gel and the light source) we observe in- 
creased light intensity in the vicinity of the crack. When 
the focal plane is slightly above the gel plane (between 
the gel and the camera), the same region becomes black. 
Hence the near tip region behaves as a diverging lens. 

Let us now consider the displacement field n(r, 6) 
around the crack tip at other 6>'s. The best defined di- 
rection is along the axis of propagation directly ahead of 
the crack (^ = 0). In Fig. 2 we present an example of the 
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FIG. 3: (Color online) (a) The stress intensity factor, K/, 
is extracted from crack opening displacement (CTOD) using 
Eq. ([!]) (circles) and from Ux (r, ^ = 0) measurements using 
Eq. (|2]) (squares). The 20% difference between the two is 
well beyond measurement uncertainties. Dashed lines are to 
guide the eye. (b) The fracture energy, F, calculated using 
Ki (derived from the CTOD using Eqs. ^ and {|3))). The 
weak increase of F with v relative to S implies that S is due 
to nonlinear elasticity, rather than dissipative processes. 



measured displacements along this symmetry axis, ob- 
tained using the particle tracking method explained ear- 
lier. We first notice that the displacements are negative, 
since particles ahead of the crack are pulled towards it. 
Comparing these measurements to the predicted LEFM 
functional form [l] 
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we obtain excellent agreement at all measured scales, 
from distances well beyond 5{v) (7mm in Fig. 2) to the 
closest near-tip vicinity that we are able to measure [9[. 
Eq. ([2]) contains the contribution of the T-stress as well, 
since the K-field alone is insufficient to describe the data. 
Thus in the region of study we have K-T dominance in- 
stead of the typically assumed i^-dominance. The next 
order term in the Ux expansion (r^/^) is the first to in- 
corporate dynamic corrections due to non-steady state 
contributions like v and Kj pLj] • We found no evidence in 
our data for the existence of such contributions. 

We are now in a position to compare the values of 
Kj extracted from the two different methods: the dis- 
placement field along the symmetry axis (using Eq. ([2])) 
and the CTOD measurements (using Eq. ([1])). We ex- 
pect to find quantitative agreement between the two mea- 
surements, as LEFM predicts that Kj is unique and has 
no angular dependence. Curves of Ki{v) extracted from 
both measurements are presented in Fig. 3(a). Although 
both describe a similar functional profile: Kj increases 
slowly until reaching a peak at v^O.bCs before decreasing 
again, surprisingly, the two Ki{v) curves quantitatively 
differ. The systematic 20% difference in the values of Kj 
is much larger than our measurement error and does not 
appear to be ^-dependent. As our measured values of 
T/E are less than a few percent, any uncertainty in their 
value could not explain this large deviation in Kj. 

Unlike the quantitative discrepancy, the qualitative be- 
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FIG. 4: (Color online) The measured strain, eyy{r,0 = 0), 
(dots) is compared to the theoretical (LEFM) prediction 
(solid curve): [gy{0 = 0,v)-Ki/E]/^ - 0.5T/E where Ki 
and T are taken from the Ux fit (cf. Fig. 2) and gy is a known 
universal function [1]. The discrepancy between the two in- 
creases with the crack velocity: (a) v = 0.20cs, (b) v = 0.53cs, 
(c) v = 0.78cs. For the higher velocity (c), LEFM predicts a 
negative strain (compression) ahead of the crack tip. 



havior of Kj is easily understood. According to LEFM, 
Kj factorizes into a universal dynamic component, which 
depends solely on the velocity and a geometric compo- 
nent, which depends on the loading and dimensions of 
the stressed object. While the first is a decreasing func- 
tion of v^ the latter increases with v through the loading. 
The non-monotonicity of Ki{v) is due to the competition 
between these two components. 

Using Ki(v)y we can now compute the material's frac- 
ture energy, F [l| 



Tiy) =K]{v)Ai{v)/E 
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where Aj{v) is a known [1] function of v. Using the values 
of Ki{v) obtained from the CTOD, we observe [Fig. 3(b)] 
that F is a slowly monotonically increasing function of 
V whose value is approximately 20J/m^. Similar values 
have been reported for other gels [l0[. F is a measure of 
the dissipation involved in crack propagation. It cannot, 
in general, be estimated from first principles, but is a 
material property which must be measured. 

We are now in a position to understand the nature 
of the processes that govern the deviations S from the 
parabolic CTOD. In principle, 5 can either be the scale of 
dissipative processes (i.e. the process zone) or it may be 
associated with nonlinear elastic processes. If the former 
were correct, one would expect r{v) and 6{v) to have 
a similar functional dependence. Comparison over the 
velocity range of 0.4 — O.Scg reveals only a modest ^30% 
increase in F, while S increases by over ~ 400%. This 
large contrast indicates that nonlinear (non-dissipative) 
elastic processes are dominant in the near-tip region at 
scales that are significantly larger than dissipative ones. 

The large discrepancy between the two values of Kj 
derived from different components of the same displace- 
ment field, motivated us to examine also the deformation 
in the y-direction at ^ = 0. Since for Mode I symmetry, 
Uy{r,0) = 0, we instead use the strain £yy{r^O = 0) = 
dyUy{r,0). In Fig. 4 we compare £yy{r^O) at three dif- 
ferent velocities, from 20% to -80% of c^, to the LEFM 



predictions [ll|. At all velocities we measure strains of 
^ 0.2-0.3, before encountering measurement limitations 
[§] when too close (~ 200/im) to the tip. LEFM pre- 
dictions for Syy{r, 0) were calculated using the Kj and T 
stress values obtained by fitting the Ux{r^ 0) displacement 
components at the same velocities, cf. Fig. 2. 

In stark contrast to the excellent fits obtained for Ux^ 
the predicted values of Syy deviate significantly from the 
measured ones. At the lower velocities, the measured 
strain increases notably more rapidly than the r~^/^ de- 
pendence predicted by LEFM [Fig. 4(a-b)], and attempts 
to fit the strains using LEFM fail. At higher veloc- 
ities [Fig. 4(c)] the deviations between measured and 
expected values are even more dramatic. 

For V > 0.73cs (for incompressible materials) LEFM 
predicts a monotonically decreasing strain which reaches 
negative values (compressive strains) as one approaches 
the crack tip. On the other hand, our measurements 
show the strain to be a positive, monotonically increas- 
ing function that reaches near-tip values that are similar 
to those measured at lower velocities. The measurements 
make more intuitive sense than the LEFM predictions - 
tensile fracture occurs under extension. 

The experiments described here present high resolu- 
tion measurements of the displacement field surrounding 
dynamic cracks. We show that the canonic theory of frac- 
ture, LEFM, fails to provide a consistent description of 
the experimental data, apparently as a result of elastic 
nonlinear it ies ^]. This, by itself, is neither surprising 
nor contradicts LEFM, as these effects may be concep- 
tually incorporated into one of its key assumptions, the 
small-scale yielding condition [1]. 

On the other hand, as fracture occurs precisely at the 
smallest scales, a description of this near-tip region is 
important. Above, we presented detailed measurements 
within the nonlinear region in this type of material. In 
addition, we demonstrated that at high v, LEFM quali- 
tatively fails to describe the tensile extensions measured 
ahead of the crack tip cf. Fig. 4(c). This underlines the 
necessity for a more complete theoretical understanding 
of essential nonlinear effects. 

Are these observations unique to elastomers or are they 
more generally valid? In the accompanying Letter [8] we 
develop a weakly nonlinear theory of the dynamic frac- 
ture of a single crack that resolves all of the discrepancies 
discussed above. Such nonlinearities are as universal as 
linear elasticity and must be experienced by any material 
undergoing fracture. Therefore, we expect these results 
to be generally applicable to any brittle material. 

The key features of brittle fracture in the gels consid- 
ered here are identical to those observed in other brittle 
amorphous materials, like glassy polymers or structural 
glasses [6j. This work, however, indicates that LEFM is 



unable to describe these features in gels, since the near- 
tip fields are nonlinearly elastic at scales encompassing 
the origin of these effects. Moreover, the near-tip dissi- 
pative processes in gels, brittle plastics and glasses are as 
different as their micro-structure. This leads us to con- 
clude that the nonlinear elastic region that bridges the 
gap between LEFM and the process zone must play a 
critical role in governing the fracture process. Thus, un- 
derstanding the dynamics within this region may be the 
key to unlocking a plethora of open questions that are 
related to the breakdown of single straight cracks. These 
include questions of stability [2| (e.g. micro-branching 
and crack oscillations), crack path selection [2], and 3D 
nonlinear focussing (e.g. crack front inertia, front wave 
nonlinear structure, and the formation of directed chains 
of micro-branches p ) . 
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